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Additive-State-Decomposition-Based Dynamic Inversion Stabilized Control for a Class 
of Uncertain Systems and Its Application to Missile Longitudinal Tracking 

Quan Quan and Kai-Yuan Cai 



Abstract 

This paper presents a novel control, namely additive-state-decomposition-based dynamic inver- 
sion stabilized control, that is used to stabilize a class of systems subject to nonparametric time- 
rs! 1 varying uncertainties with respect to both state and input. By the additive state decomposition and 
a new definition of output, the uncertain system is transformed into an uncertainty-free system 
with relative degree one, in which all uncertainties are lumped into a new disturbance at output. 
Subsequently, dynamic inversion control is applied to reject the lumped disturbance. Performance 
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analysis of the resulting closed-loop dynamics is given that the stability can be ensured if a key 
parameter is sufficiently small. Finally, this proposed control is applied to a tracking problem for the 
missile's angle of attack. Numerical simulation demonstrates its effectiveness. 

I. Introduction 



The stabilization problem for control systems subject to uncertainties with respect to both state 
and input is a challenging problem, which has become a subject of keen interest in recent years. It is 
well known that uncertainties at input are often present in practical systems, for example input with 
an uncertain gain, or nonlinearities inherently aroused from dead zone, quantization and backlash, 
etc. The existence of uncertain inputs is a source of degradation or instability in the performance of 
the system. In addition to uncertainties in the input, there exist plant uncertainties, which originate 
from various sources, such as variation of plant parameters, inaccuracy arising from identification, 
etc. Since aerodynamic parameters are functions of the flight condition, a missile is a nonlinear and 
rapidly parameter-varying plant due to operation in a large range of aerodynamic conditions. The 
aerodynamic parameters exist in state matrix and input matrix _Q|-|__, and are uncertain in practice. 
According to this, missile control problem can be formulated into a problem of stability of robust 
control systems subject to uncertainties with respect to both state and input. 

In this paper, the stabilization control problem for a class of systems subject to nonparametric time- 
varying uncertainties with respect to both state and input is investigated. Before introducing our main 
idea, some accepted control methods in the literature to handle the uncertainties are briefly reviewed. 
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A direct way is to estimate all of the unknown parameters, then use them simultaneously. Lyapunov 
methods are used to analyze the stability of closed-loop system. In Q, L\ adaptive control architecture 
is proposed for systems in the presence of unknown input gain, time-varying unknown parameters and 
disturbances. This controller leads to uniform performance bounds for the system's input and output 
signals, which can be systematically improved by increasing the adaptation rate. In |4), the problem 
of a planar rigid body, with unknown rotational inertia and an unknown input nonlinearity, tracking a 
desired angular velocity trajectory is addressed using adaptive feedback control. By parameterizing the 
input nonlinearity, an adaptive feedback linearization controller is developed. In 0, two asymptotic 
tracking controllers are designed for output tracking of an aircraft system in the presence of parametric 
uncertainty and unknown, nonlinear disturbances, which are not linearly parameterizable. An adaptive 
extension is then presented, in which feedforward adaptive estimates of the input uncertainty are used. 
In Q-[5], an adaptive controller may require many integrators corresponding to many unknown 
parameters for an uncertain system because each unknown parameter needs an integrator to estimate. 
This will lead to a resulting closed-loop system with a reduced stability margin. In addition, the 
estimates may not approach the true parameters without persistently exciting signals, which are 
difficult to generate in practice especially when the number of unknown parameters is large |6] 
pp. 111-118]. The second way is to design an inverse control by a neural network to cancel the 
input nonlinearities in order to obtain a linear function [7]], [HI. Contrary to the traditional inverse 
control schemes, a neural network is employed to approximate an unknown nonlinear term instead 
in [9]. In fact, the neural network methods can also be considered as an adaptive control method just 
with different basis functions. So, this way also has the same problems as that for adaptive control. 
The third way is to adopt sliding mode control due to its inherent advantages of fast response and 
insensitivity to plant parameter variation and/or external perturbation. In ifTOll . a new sliding mode 
control law based on the measurability of all of the system states is presented to ensure the global 
reaching condition of the sliding mode for systems subject to uncertainties with respect to both state 
and input. Based on it, an output feedback controller was further proposed in ifTfl . Whereas, sliding 
mode controllers essentially rely on infinite gains to achieve good tracking performance, which is not 
always feasible in practice. The drawbacks of the high-gain feedback solutions are related to the fact 
that they may saturate the joint actuators and excite high-frequency modes, etc. 

In order to overcome these drawbacks, this paper proposes a novel control, namely additive-state- 
decomposition-based dynamic inversion stabilized control, which is used to stabilize a class of systems 
subject to nonparametric time-varying uncertainties with respect to both state and input. Unlike many 
literatures, uncertainties are nonparametric and time-varying. Therefore, the existing methods, namely 
adaptive control method and neural network method, are not applicable directly. The key idea behind 
the proposed method is to lump nonparametric time-varying uncertainties with respect to both state 
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and input into one disturbance by additive state decomposition, and then to compensate for it. The 
proposed controller is continuous and can produce asymptotic stability if the uncertainties are time 
invariant. At this point, it overcomes the drawbacks of sliding mode controllers. Finally, the additive- 
state-decomposition-based dynamic inversion stabilized controller is applied to tracking control for 
the missile's angle of attack. Unlike the control methods proposed in [l]-[2j, less assumptions are 
imposed on the time-varying term. Furthermore, the design of the controller is easier and the structure 
is simpler. 

II. Problem Formulation and Additive State Decomposition 

A. Problem Formulation 

Consider a class of systems subject to nonparametric time-varying uncertainties with respect to 
both state and input as follows: 

x (t) = A x (t) + b(h (t, u)+a (t, x)) ,x(0) = x (1) 

where x (t) G R n is the system state (measurable), «(t)el is the control, b G R n is known constant 
vector, Aq G R nxn is a known matrix, h : [0, oo) x R — > R is an unknown nonlinear function, and 
a : [0, oo) x R n — > R is an unknown nonlinear time-varying disturbance. For system (Q]), the following 
assumptions are made. 

Assumption 1. The pair (Ao,b) is controllable. 

Assumption 2. The unknown nonlinear function h satisfies h(t, 0) = 0, 1^1 < \u\ , || > 

and | HI < l hu , Vn G R, Vi > 0, where l ht ,L hu Jh u > 0. 

Assumption 3. The time-varying disturbance a (t, x) satisfies \a (t, x)\ < k a \\x (t)\\+5 <r (t) , ||ff|| < 
l ffm and |f | <l at Wx^W+d^t), where % = [^--- |]eK lxn , K,6 a (t) J^l^d* (t) > 0, 
Vt > 0, and 5 a (t) , (t) are bounded. 

{u \u\ > \i 
, then it can be 
\u\ < \i 

reformulated as h(t,u) = u + 8(t) , where \5 (t)\ < fj,. In practice, we do not need to know the 
parameters l ht , l_ hu , l hu ,k a ,5 a ,l r7x ,l C7t ,d C7 . 

The control objective is to design a stabilized controller to drive the system state such that x (t) — > 
as t — > oo or it is ultimately bounded by a small value. In the following, for convenience, we will 
drop the notation t except when necessary for clarity. 

B. Additive State Decomposition 

In order to make the paper self-contained, additive state decomposition in jT2l. lfT3l is recalled 
briefly here. Consider the following 'original' system: 

f(t,x,x) = 0,x(0) = x (2) 
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where x S W 1 . We first bring in a 'primary' system having the same dimension as (0, according to: 

fp(t,Xp,x p ) = 0, x p (0) = x pfl (3) 

where x v € ffi". From the original system (0 and the primary system ([3]) we derive the following 
'secondary' system: 

/ (*, x, x) - f p (t, x p , Xp) = 0, x (0) = x (4) 
where x p € M n is given by the primary system Q. Define a new variable x s € M n as follows: 

x s — x Xp. (5) 
Then the secondary system © can be further written as follows: 

/ (t, x s + x p , x s + x p ) - f p (t, x p , x p ) = 0, x s (0) = x - x P)0 . (6) 
From the definition ([5]), we have 

x (t) = Xp (t) + x s (t),t> 0. (7) 

Remark 2. By the additive state decomposition, the system © is decomposed into two subsystems 
with the same dimension as the original system. In this sense our decomposition is "additive". In 
addition, this decomposition is with respect to state. So, we call it "additive state decomposition". 

As a special case of ©, a class of differential dynamic systems is considered as follows: 

x = f(t,x),x(0) = sc , 

y = g(t,x) (8) 

where and y G R m . Two systems, denoted by the primary system and (derived) secondary 

system respectively, are defined as follows: 

Xp — fp (t, Xp) , Xp (0) — Xpfl 

V v = 9 P {t, x p ) (9) 

and 

X s — f (t, Xp -\- X s ) fp (t, Xp) , X s (0) — Xq Xpfi, 

y s = g{t,x v + x s ) - g P (t,x p ) (10) 

where x s = x — x p and y s = y — y p . The secondary system (fTOl) is determined by the original system 
© and the primary system (O. From the definition, we have 

x (t) = xp (t) + x s (t) ,y(t) = y p (t) + y s {t),t> 0. (1 1) 
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III. Additive-State-Decomposition-Based Dynamic Inversion Stabilized Control 

In this section, the considered uncertain system is first transformed into a determinate system but 
subject to a lumped disturbance by additive state decomposition. Then a dynamic inversion method 
is applied to this transformed system. Finally, the performance of the resulting closed-loop system is 
analyzed. 

A. Decomposition 

Since the pairs (Aq, b) is controllable by Assumption 1, we can always find a vector k G W 1 such 
that A + bk T is stable. This also implies that there exist < P,M G R nxn such that 

P(A + bk T ) + (A + bk T ) T P=-M. (12) 

According to this, we rewrite the system £T|) to be 

x = Ax + b (h (t, u) - k T x + a (t, x)) , x (0) = x (13) 

where A = A$ + bk T . Consider the system (fl"3l as the original system. We choose the primary system 
as follows: 

x p = Ax p + bu, x p (0) = 0. (14) 

Then the secondary system is determined by the original system ([T3l and the primary system (fl4l 
with the rule (TTOt . and we obtain that 

x s = Ax s + b (—u + h (t, u) - k T (x p + x s ) + a (t, x p + x s )) , x s (0) = x - (15) 

According to (fTTl) . we have 

x = x p + x s . (16) 

The secondary system is further written as 

x s = Ax s + b (—u + h (t, u) — k T x + a (t, x)) , x s (0) = Xq. (17) 

Before proceeding further with the development of this paper, the following preliminary result is 
needed. 

Theorem 1. Suppose det(sJ„ — A) = (as + 1) (/3 n _is n_1 + • • • + P\s + /3o) • If Assumption 1 
holds and (sl n — A)^b is asymptotically stable, then there exists a c G M n such that 

1 



c T (s/ n - A) _1 6 



as + 1 



where a > 0, c = (A r " 1 ) T c G R n , c = $ n _ x ... fa fa 
defined that adj(s/ n - A)b = N^- 1 ■■■ 1} T . 
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Proof. Given a vector c € R n , we have 



c T (si n -A)^b= c y^-^ b 

y ' det(sJ n - A) 

c T N[s n ~ l 



(as + 1) {Pn-lS n - 1 +... + p xS + fay 

If the pair (Aq, b) is controllable by Assumption 1, then the pair (A, b) is controllable as well. Conse 
quently, the matrix N is nonsingular 1 14, Lemma 3]. Let c = (iV -1 ) T c€lR n ,c = 
R n . Then 



J ^— ro " - [ /?„_! ••• /3i A, fe 



r/ r A»-l* n_1 + • • • + Pis + A) 



(as + 1) (/^-is™- 1 + ■ • • + fas + fa) " 
Since (s/ n — A)~ 1 b is asymptotically stable, f3 n -is n ~ 1 + • • • + f}\s + /3o is an asymptotically stable 
(n — l)th order polynomial. So, we have c T (sI n — A)~ x b = a ^ +l . □ 

Theorem 1 shows that we can always find aceR" such that G = a ^ +1 . According to Theorem 
1, we define a new output y = c T x and rearrange (fT4b - ([T6l) to be 

x p = Ax p + bu, x p (0) = 
y = c T x p + di (18) 

where di = c T x s = G (—u + h (t, u) — k T x + a (t, x)) + c T e At xo is called the lumped disturbance. 
Furthermore, (fT8l ) is written as 

y = Gu + d l . (19) 

The lumped disturbance d\ includes uncertainties, disturbance and input. Fortunately, since y p = 
c T x p = Gu and the output y are known, the lumped disturbance di can be observed by 

di = y - Gu. 

It is easy to see that d\ = di. 

B. Dynamic Inversion Control 

So far, by the additive state decomposition, the uncertain system (Q]) has been transformed into a 
determinate system (fT9l but subject to a lumped disturbance, which is shown in Fig.l. 

For the system (fT9l . since G is minimum phase and known, the dynamic inversion tracking 
controller design is represented as follows: 

u = -G- l d h (20) 

Substituting (f20]> into (O results in 

y = -GG- 1 ^ + dl 

= -di + di = 
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P:x=A x+b(-) G = c T (sI n -A)~ l b 




(b) Determinate system subject to a lumped disturbance 



Fig. 1. Model Transformation 

where d = di is utilized. As a result, the perfect tracking is achieved. However, the proposed controller 
(l2Qb cannot be realizable. By introducing a low-pass filter Q, the controller (1201 ) is modified as follows: 

u = -QG- l d h (21) 

By employing the controller (I2TI) . the output becomes 

2/=(l-Q)4 (22) 

Since Q is a low-pass filter and the low-frequency range is often dominant in a signal, it is expected 
that the output will be attenuated by the transfer function 1 — Q. However, the mention above is only 
a straightforward explanation. A detailed analysis is given in the following section. 

C. Performance Analysis 

Since the lumped disturbance di involves the input u, the resulting closed-loop system may be 
instable. Next, some conditions are given to ensure that the control input u is bounded. Substituting 
d\ into (1211) results into 

u = -QG' 1 [G (-u + h (t, u) -k T x + a (t, x)) + c T e At x ] 



= Q(u-h{t,u) + k T x - a (t, x)) - QG- 1 c T e At x . (23) 
Multiplying Q^ 1 on both sides of d23l yields 

Q' x u = u - h (t, u) + k T x - a (t, x) + ^ (24) 
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where £ = — G~ l c r e At XQ. Since G = a ^ +1 and A is stable, the term £ will tend to zero exponentially. 
A simple way is to choose Q = jjxj, where e > can be considered as a singular perturbation 
parameter. By the filter Q, (l24l) is further written as 

eii = -h (t, u) + k T x - a (t, x) + £. (25) 

So far, the dynamics (fT3T > and (1251 ) form a new closed-loop dynamics. In geometric singular per- 
turbation theory Q3), the behavior of singularly perturbed systems is determined using geometric 
constructs of the reduced-order models, which are obtained by substituting e = in (l25l ). According 
to this, we have 

« -h (t, u) + k T x - a {t, x) + £. 

By this result, (IT3T > becomes 

x ?» Ax + &£, x (0) = xq. 

The dynamics above imply that x — >■ as t — > oo if A is stable. Furthermore, it is bounded by (|25T ). 
The following theorem will give an explicit bound on e which can ensure the stability of closed-loop 
dynamics forming by (fT3T ) and (|25T ). 

Theorem 2. Suppose i) Assumptions 1-3 hold, ii) the controller is designed as (1211 ) with Q = ^W, 
iii) e satisfies 

6 < 27 s T72 (26) 



where 



70 = A min (M) 

71 = 2(^11+^)11611 + 2^ 

72 = ll^ll ||b|| + \\A\\ (\\k\\ + l a J + + k a . (27) 



Then the state of system £[]) is uniformly ultimately bounded with respect to the bound y (f 11 "^' 7 + 

Furthermore, if ^ t 5(j — > and d CT — > 0, then the state x — > as t — > oo. In addition, if e is sufficiently 
small, then the uniformly ultimate bound of state x is y / 2X "™*( P A fj^^o- + <^o-J • 
Proof. See Appendix. 

Remark 3. According to Theorem 2, if the uncertainties are time invariant, namely ^ t = and 
da = 0, then the proposed controller can produce asymptotic stability. Also, from Theorem 2, a 
sufficiently small e will satisfy d26l) and the ultimate bound will be reduced by decreasing e. However, 
we should point out that a small e in turn will result in a reduced robustness of the closed-loop system. 
For example, consider a simple situation where h (t, u) = u(t — r) . Given a r > 0, the dynamical 
system eu = — u (t — t) will loss stability by choosing a sufficiently small e no matter how small the 
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delay r Therefore, an appropriate e > should be chosen to achieve a tradeoff between tracking 
performance and robustness. We do not need to know the parameters 1^,1^,1^,^, 5 c ,l c<ii ,l Ct ,d a 
and choose e according to practical condition. 

From the above analysis, the design procedure is summarized as follows. 

1. Design a state feedback gain feel" such that A = Aq + bk T is stable with det(sl n — A) = 
+ /3is + /3 ) . 

1] T and c 



(as + 1) (Pn~is n - 1 + 



2. Design an output matrix c 

i T 



[N- l ) T c € R n where adj(sl n - A)b = N^- 1 • • • 1lT 



i 



ft A) 



3. Design d2B with G = -4rr and Q = -4r- 

4. Choose an appropriate e > in practice. 

IV. An Application to Missile Longitudinal Tracking 

The following linear parameter-varying system approximates a missile's longitudinal axis dynamics 
ID-El: 

-Z a (t) 1 
-M a (t) 



a (t) 




m _ 







a(t) 









+ 










_ M S (t) _ 



(S m (t) + d(t)) 



(28) 



where a is the angle of attack, q is the pitch rate, and S m is the fin deflection. We assume that a, q 
are measurable and Z a (t) is knowrj^. In the following we will formulate (|28T ) as dTJ. The objective 
is to let a (t) — > r (t) as t — > oo, where r (t) is known and twice continuously differentiable. 



A. Problem Formulation 
Define 

We have 



xi 

X 2 



,x\ = a — r,X2 



-Z a a + 



xi = x 2 

X2 = —Z a a — Z a a + q — r 

= -Z a a - Z a {-Z a a + q) — M a a - r + Ms (5 m + d) 



(29) 



'The characteristic equation of eii = —it (t — r) is es + e~ ST = 0, which can be approximated by (e — r) s + 1 = 0. 
Therefore, one solution of the characteristic equation is s ~ 737- If £ < t, then s > 0, namely the dynamical system 
eii = —u (t — t) is unstable. 

2 The requirement of an exact Z a (t) is only used for the problem formulation. In practice, the parameter can be inexact. 
This will be shown in the following simulation. 
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In practice, we only know the approximation of the parameters M a (t) , Ms (t), namely M a (t) , Ms (t) 
With these information, the control 5 m is designed as 



5 m = t>- (k T x + Z a a + Z a ( — Z a a + q) + M a a + r + u 
M & V 



(30) 



where the additional part k T x is to stabilized the closed-loop system, while the other additional part 
[z a a + Z a (—Z a a + q) + M a a + r^j is to compensate for — Z a a — Z a (—Z a a + q) — M a a — r. 
By substituting (1301) into d29l) , the missile's longitudinal axis dynamics are formulated as ([T]), namely 






1 











x + 












1 



(/»(*,«) + <r(t,x)) 



(31) 



where h (t, u) = j^u and a (t, x) = ~± (k T x + Z a a + Z a {-Z a a + q) + M a a + rj - Z a a - 
Z a (—Z a a + q) — M a a —f +Msd. Obviously, if M a = M a and Ms = Ms, then the system above 
becomes 

oil r o 

x = x+ (it + M s d) . 

ki k 2 1 

B. Controller Design 

For the system (l3Tb . according to the design procedure at the end of Section III.C, we have the 
following design. 



Step 1. According to the design procedure, we design k 

1 



-12 -7 



resulting in 



A = 

with det(s/ 2 ~A) = (\s + l) (4s + 12) . 

Step 2. According to A and b, we can obtain c 





-12 



-7 



4 12 



and N 



the output matrix is designed as c 



12 4 



T 



1 

1 



Consequently, 



Step 3. Design controller it = — QG 1 di, where G = T^pj, Q = ^pjj = c T x — Git and e will 
be specified latter. 

C. Numerical Simulation 

Similar to [2], we choose the parameters 

Z a (t) = 2.35 - 2.05 sin(0.05vrt),M a (t) = 5 + 375 sin(-0.014vrt) 

M s (t) = 3-2.5 cos(0.04vrt), d(t) = 0.1 (32) 

where the working interval < t < 100s. The initial conditions are a (0) = — 1 and q (0) = 0. Denote 

-Z n (t) 1 



F 



-M a (t) 



. When t = 0s, eigenvalues of matrix F are — 1.1750±1.9025i. When t = 50s, 
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eigenvalues of matrix F are ±17.4244. The considered linear parameter-varying system is unstable 
during the working interval. In practice, we assume that we only know the approximation of the 
parameters, namely, 

Z a (t) = 2.2 - 2sin(0.057rf),i Q (t) » 

M a (t) = 4 + 360sm(-0.0147rf),M 5 (i) = 2.8 - 2.1 cos(0.04vrf) (33) 

where Z a (t) is approximated to Z a (t). 

First, according to the known information, we design 8 m but without u, namely 

^ m = ~W ^ x + + ^ a (~^ a01 + q ) + ^ aOL + ^) • 

The simulation is shown in Fig.2. The divergence suggests that feedback stabilization and compensa- 
tion with approximated information may not achieve the tracking task, especially when the parameters 
vary in a large range. Let us add the control term u into S m with e = 0.06. The simulation is shown 
in Fig. 3. As shown, the controller with the stabilized control term can make a track the reference r 
with a small tracking error in the presence of time-varying uncertainties. This is consistent with the 
conclusions of Theorem 2. 




-2' 1 1 1 1 1 1 1 1 1 1 

2 4 6 8 10 12 14 16 18 20 



.1 

!i 
i 

i ■ 


... r 

a 

- q 







10 20 30 40 50 60 70 80 90 100 

time(sec) 

Fig. 2. Tracking performance without additive-state-decomposition-based stabilized control term 

In order to show the asymptotic convergence, we assume that the tracking reference r (t) = 1 and 
these time-varying parameters are frozen at t = 50s, namely Z a (t) = 0.3, M a (t) = —298.4, Mg(t) = 
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1.5 




10 20 30 40 50 60 70 80 90 100 

time(sec) 



Fig. 3. Tracking performance with additive-state-decomposition-based stabilized control term in the presence of time- 
varying uncertain parameters. 

0.5, d(t) = 0.1. Similarly, we assume that we only know the approximation of the parameters 
Z a (t) = 0.3, Z a (t) « 0,M a (t) = -287, M$ (t) = 0.7. By the same controller with additive-state- 
decomposition-based stabilized control term with e = 0.03, we get the simulation result shown in 
Fig. 4. The tracking error tends to zero. This result is consistent with the conclusions of Theorem 2. 

V. Conclusions 

Stabilized control for a class of systems subject to nonparametric time-varying uncertainties with 
respect to both state and input is considered. Our main contributions lie in two transformations: i) the 
presentation of a new decomposition scheme, called additive state decomposition, which transforms 
the uncertain system into an uncertainty-free system; ii) and a definition of output matrix which 
transforms the uncertainty-free system into a first order system. Based on these transformations, 
nonparametric time-varying uncertainties with respect to both state and input is lumped into a new 
disturbance at the output of the first order system. This can explain why the proposed method can 
handle these nonparametric time-varying uncertainties. The proposed control scheme has two salient 
features: less information of one system required and simple design procedure with less tuning 
parameters. 
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2.5 



3.5 



4.5 




Fig. 4. Tracking performance with additive-state-decomposition-based stabilized control term in the presence of constant 
uncertain parameters. 



VI. Appendix: Proof of Theorem 2 
Denote v = h (t,u) — k T x + a (t, x) . Then the system (TT3l becomes 

x = Ax + bv 



and the derivative of ev is presented as 



dh dh 

ev= d- u eu + e m 



. rp , da . da / —— . 
ek x + —ex + e— ( (1251) is used) 

ox at 



dh 
du 



(-'' + () + '( -A: T + |^ J (Ac + bv) + < ~ + r 



9/i 

at 



da 
dt 



dh 



+ e I -k T + ] b 



da 



du \ dx 
Consequently, the closed-loop dynamics (1251 and (fTBl are 



5/i 9<T e?/i ^ / , t da., 

v + e m +e m + ^ + e [- k + d~ x ^ Ax - 



x = Ax + 6v 



1 dh f T da 
e du \ dx 



^ dh ^ da ^ 1 dh ^ / _|_ ^°"\ j 



(34) 



Design a candidate Lyapunov function as follows: 



V = x T Px + v 2 
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where < P £ M nxn satisfies (fl2l ). Taking the derivative of V along the solution of (1341) yields 

V" = x T (PA + A T P) x + 2x T Pbv + 2v 



ldh f ^ T da . ^ 
e 3w V cte 



J. 9<t\ , ^ f dh da ldh 



t / a at n ( dh da 1 dh A 

I T {PA + A T p)x + 2v ^_ + _ + __ t j 



ldh ( ^ T da . ^ 
e du \ dx 



v 2 + 2x T 



Pb -A T (-k T + ^-\ 



dx J 



v. 



By (O, we have x T (PA + A T P) x < -70 ||x|| 2 , where 70 = A min (M) . Then 
V < -70 IN| 2 + 2 |v| |u| + l at \\x\\ + d a + ^ |f| 

2 ( — —71 J + 272 |f I ||x|| (Assumptions 2-3 are used) 



(35) 



Next, before proceeding further, we will derive the relationship between u and u. By the mean-value 
theorem, we have 



h(t,u) = h(t,0) 



dh{s) 



ds 



where < c < u. Furthermore, since v = h(t,u) — k T x + a (t, x), we have 



/ dh (s) 
U= { — 



-1 



(v + k T x - a (t, x) -h(t,0)) . 



Further by Assumptions 2-3, the equation above becomes 



M < — {\v\ + 
<t^(H + 



x\\ + \a (t, x)\) 

+ ker) \\x\\ + 6 a 



With this inequality above and relationship 

2 | V | ( l -^s a + da + l J±\z\)< \y + J- ( l J±5 a + d a + ^ lei 



the inequality (1351 ) becomes 



V < -70 INI 2 - ( \l hu - 7i) v 2 + + d a + ^ |e| ) + 2 ( 72 + Z fft ) [ v | . 

2 , 1,2 



Since the inequality 2ab < a + b always holds, then 



2 (72 + W H INI < ^\\A\ 2 + -{i2 + i* t ?v 2 . 

2 70 



Consequently, we have 



T > ^ 70 || i|2 . ^ ( ht c . , ^„ it 1 

^<-— ||x|| +— — ^ + + — Id 



-lh u -11- — (72 +^J 2 

e 7o 



v 2 . 
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Furthermore, we have 

V<- V (e)V + ^-(^5 a + d a + l -f\t\ 

where rj (e) = min( 2A 70 ( P ) , \lh n ~ 7i — ^ (72 + ^o-J 2 )- If <l26l) is satisfied, then 77 (e) > 0. So, we 
have F (i) -> B (^yif- (7^ + <k + as t -4 00, namely 

l|i(t>l| -K&&(b +d » + T^)) 

as t -> 00, where = {£ E < 5}, 5 > 0. The notation x(t) -4 means min 

2/e8(<5) 

- y| -4 as t -4 00. Since £ -)• as t ->• 00, ||x(t)|| ->■ B O Amin (p)^( ^~ (iT^ + • 
Furthermore, if lh t S a — > and d a — > 0, then the state x — > as i — > 00. In addition, if e is sufficiently 
small, then 77 (e) = 2A 70 ^ , namely the uniformly ultimate bound of state x is 2Xmax \ p ) -±- ( ^-<5 CT + d a 



— — (7^) ' llallltl J Lllt uiiiiuiim^ uiuiiiaLt uuunu ui aiaLt j> is — j 
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